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Abstract 

The simulation of the interaction of two simultaneously propagating air streamers 
of the same polarity is presented. A parametric study has been carried out using an 
accurate numerical method which ensures a time-space error control of the solution. For 
initial separation of both streamers smaller or comparable to the longest characteristic 
absorption length of photoionization in air, we have found that the streamers tend to 
merge at the moment when the ratio between their characteristic width and their mutual 
distance reaches a value of about 0.35 for positive streamers, and 0.4 for negative ones. 
Moreover it is demonstrated that these ratios are practically independent of the applied 
electric field, the initial seed configuration, and the pressure. 

Electrical breakdown in air gaps often involves the development of fast ionizing waves 
that take the form of thin filaments called streamers [1] . These complex and highly nonlin- 
ear phenomena precede the formation of sparks and leaders. Similar filamentary structures 
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are observed in sprites: large-scale discharges which appear at low pressure in the altitude 
range of about 40 to 90 km above large thunderstorms 2, 3 . Streamers often form branches. 
Then the individual streamer heads carry charges of the same polarity, and thus they are 
electrostatically repelling each other. However, careful observations show also the opposite 
behavior: streamers may even merge or reconnect [H[51[B]. So far, only few numerical studies 
have been carried out on streamer interactions. First, the impact of charges of numerous 
positive streamers, propagating simultaneously, on the electric field and on the velocity of 
one single positive streamer was studied in [7] for a 2D configuration. Later on it was shown 
in [8] that the electrodynamics of a 2D array of negative streamers substantially differs from 
the one of a single streamer due to their electrostatic interaction. Recently, the interaction of 
two streamer discharges in air and other oxygen-nitrogen mixtures has been studied for a 3D 
cylindrical configuration in [9] . These authors showed that two competing mechanisms have 
to be considered: attraction due to nonlocal photoionization between streamers and electro- 
static repulsion of streamer heads due to their space charges of same polarity. Qualitatively, 
they concluded that streamer merging is favored if the distance between streamers is smaller 
or comparable to the longest absorption length of photoionization. However, a more precise 
or even quantitative description of streamer merging remains an open problem. 

In this work we carry out a parametric analysis of the dynamics of interacting streamers 
to better understand the conditions of streamer merging. Therefore we simulate the inter- 
action of two streamers of the same polarity by means of an accurate numerical method 
which ensures a time-space error control of the solution. Taking into account that fully 3D 
simulations are still computationally expensive to carry out an extensive parametric analysis, 
we have considered a Cartesian 2D geometry to reduce computational costs. This study is 
then a first step towards a quantitative analysis of streamer merging in real 3D configura- 
tions. However, in agreement with [EJUS], we consider that the main characteristics of the 
dynamics of interacting streamers in two dimensions, as described below, will be similar in 
three dimensions. In practice, we aim at defining geometrical parameters independent of 
the multivariable complex physical settings and thus a simple characterization that may be 
extended to more general configurations. 

We consider the classical fluid model for air given by drift-diffusion equations self-consistently 
coupled with Poisson's equation 112 ) : 



d t n c + V • (n c v c ) - V • (D c Vn e ) = n c a\v c \ - n c r]\\- c \ - n c n p f3 cp + n n j + S ph , 
d t n p + V • (n p v p ) - V ■ (D p Vn p ) = n c a|v c | - n c n p /3 cp - n n n p /3 np + S ph , 
d t n n + V • (n n v n ) - V ■ (D n Vn n ) = n c ^|v | - n n n p /3 np - n„7, 




e V • E = -q c (n p - n n - n e ), E = -VV, (2) 

where nt is the density of charged species i (e: electrons, p: positive ions, n: negative ions), V 
and E stand, respectively, for the electric potential and field, and Vj = /i;E is the drift velocity. 
We denote by Di and fa the diffusion coefficient and the mobility of charged species i, q c is 
the absolute value of the electron charge, and £o is the permittivity of free space. Moreover 
a is the impact ionization coefficient, rj stands for the electron attachment coefficient, /3 cp 
and /3 np are, respectively, the electron-positive ion and negative-positive ion recombination 
coefficients, and 7 is the detachment coefficient. All these coefficients depend on the local 
reduced electric field E/N, and thus vary in time and space, where E = |E| is the electric 
field magnitude, and ./V is the air neutral density. Further details are given in |15j . For <S p h, 
the photoionization source term for air, we use the 3-Group SP3 model derived in [13] with 
Larsen's boundary conditions |14j . In this model, three equivalent absorption lengths are used 
to model the wavelength dependence of the photoionization source term in air. At ground 
pressure these absorption lengths are (Aipo 2 ) _1 = 0.1408cm, (A2PO2) -1 = 0.0561cm, and 
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(X3PO2) 1 = 0.0105cm, where po 2 = 158.9 Torr is the partial pressure of molecular oxygen 
in air. 

To handle the stiff modeling equations, we have recently developed in [15] a numerical 
scheme for multi-scale streamer discharge simulations for general Cartesian multi-dimensional 
geometries. The strategy is based on a second order time adaptive integration with a split- 
ting technique and dedicated solvers, and space finite volume multiresolution for dynamic 
grid adaptation. These features involve important efficiency gains in terms of CPU time and 
memory space while ensuring a time-space error control of the solution. Poisson's equation 
and Helmholtz equations issued by the photoionization model are considered on the adapted 
grid, and the resulting sparse linear systems are solved with MUMPS [161 ITT] . In all simula- 
tions carried out in this work (except for some computations in figure [5] where a finer grid is 
used) a time-space accuracy tolerance of 10~ 4 is used for a space resolution of 3.9 fim. This 
choice guarantees a sufficiently fine time-space representation of the physics, and numeri- 
cal results disclosing practically the same physical behavior for higher spatial resolutions and 
tighter accuracy tolerances. The domain size is carefully chosen for the various configurations 
and with the same space resolution of 3.9 /jm such that no interference of the boundaries is 
evidenced. 

In order to preserve similarity constraints at different pressures (corresponding to dis- 
charges at different altitudes), all variables are scaled by an appropriate power of the ra- 
tio of the air density at ground pressure N , and the air density N at the given pressure 
[T5] . For example, the space dimensions scale as x = xqNq/N, the charge densities as 
n = 72 N 2 /Nq, the electric field as E = E N/N , and the time as t — t N /N, where the 
subscript "0" refers to the values at ground pressure. Our reference (ground) density for air 
is N Q = 2.688 x 10 19 citT 3 at a temperature of 273 K. 

We study two separate discharge filaments initiated by placing two identical Gaussian 
plasma clouds into a homogeneous background electric field generated by two remote elec- 
trodes. The pressure, streamer polarity, applied electric field, and the configuration parame- 
ters of the initial Gaussian seeds (the width a, the maximum density n max , and the mutual 
separation yo) are varied in order to derive a condition for the merging of both streamers. 

Figures QJa)-(c) show snapshots of different phases of the evolution of the net charge 
density (p = n p — n c — n n ) for two simultaneously propagating positive discharge filaments 
at ground pressure (N/N n = 1.0) for an initial seed separation of yo = 0.2Nq/N cm, a 
maximum density of n max = 10 13 N 2 /Nq cm -3 , a width a — 0.02No/N cm, and an applied 
electric field E a = 1.5i?bd, where i?bd = 32 N/No kV cm~' is the breakdown electric field 
in air. The complete trajectory of the maximum electric field is also shown. Following 
this trajectory, the transition time t tT , between the streamer repulsion and the beginning of 
their merging, is defined as the instant when the separation of the two trajectories starts 
to decrease. For the condition of Figure [1] we have tt r = 6.3 ns. Figure [2ja) shows at this 
time the distribution of the net charge density. We define in Figure [2jb) the mutual distance 
of both filaments 8 by the separation of the maximum net density p max in each head; and 
the streamer head characteristic width T> as the planar length set by the cross section in the 
direction perpendicular to the applied electric field E a , where the charge density is locally 
higher than half the maximum net density: p > p max /2. 

We have found that the merging condition for two streamers in air could be represented by 
the value of the ratio V/8 which appears to be very stable for a given polarity of discharges. 
To illustrate this, Figure [3] shows the time evolution of the ratio V/8 for two positive or two 
negative streamers. In this figure, initial separation of the seeds yo = 0.15 Nq/N cm and 
0.2iVo/iVcm are considered, and other parameters are similar to Figure 1. First, Figure [3] 
shows the time evolution of the ratio T>/8 for two negative streamers as these discharges 
may propagate without photoionization source term. When photoionization is included, the 
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cm cm cm 

Figure 1: (color online) Time evolution of the net charge density for two positive streamers 
at ground pressure (N/Nq = 1.0) for an applied electric field E a = 1.5-Ebd and two Gaussian 
seeds with y = 0.2 N /N cm, n max = 10 13 N 2 /N$ cm" 3 , and a = 0.02iV /iV cm. (a) t = 
5.0 ns: well-developed streamers repulsing each other, (b) t — ttr — 6.3 ns: transition between 
repulsion and merging, (c) t = 8.0 ns: propagation of a single discharge. Black solid line: 
trajectory of the maximum electric field. 
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Figure 2: (color online) (a) Net charge density and (b) Region where the net charge density 
is higher than half the maximum net charge density at t = ttv for the same condition as in 
Figure [ljb). Definition of the streamer head characteristic width T> and the mutual distance 
of both filaments 8. 
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ratio V/5 increases and reaches a value of 0.4 at the moment of transition (4 ns in case of 
yo = 0.15 No/N cm, and 5ns in case of yo = 0.20 No/Ncm) between streamer repulsion and 
the beginning of their merging, and streamer merging is then observed. It is interesting to 
note that when photoionization is omitted, the ratio V/5 is less than 0.4, and no merging 
is evidenced. Conversely photoionization must be included for positive streamers as it is 
necessary for their propagation. In this case Figure [3] shows that the ratio V/5 increases 
and reaches a value of 0.35 at the transition time 5.6ns for yo — 0.15 No/N cm, (6.3ns for 
yo = 0.20 No/Ncm), and then streamer merging is observed. These results clearly show the 
crucial role of the photoionization source term on streamer merging. In the following we show 
that the value of the ratio V/5 at the transition time depends on the polarity of discharges 
but is practically independent on the multivariable physical settings. For instance, for the 
same parameters of the seeds as in Figure [T] (n max = 10 13 N 2 /Ng cm~ 3 , yo = 0.2 No/Ncm, 
and a = 0.02 No/N cm) we have performed a set of simulations with an applied electric field 
E a varying from 1.5i?bd to 2.5i?bd at ground pressure. Figure |4] shows that the resulting 
values of V/5 at the transition times for merging streamers are roughly constant and of 
about 0.35 for positive polarity discharges, and approximatively 0.4 for negative ones. Hence 
V/5 is independent of the applied electric field E a for both streamer polarities. Similarly, 
Figure [5] shows the dependence of V/5 on n max , the maximum density of seeds, for a fixed 
applied electric field E a — 1.5i?bd- Same previous ratios are obtained whereas T>/5 is slightly 
decreasing as rt max decreases for positive streamers. However a variation of four orders of 
magnitude for n max involves only a variation of about 15% of T>/5. In the same way we have 
tested the variation of V/5 with altitude (not shown here). We have considered altitudes 
ranging from the ground up to 80 km, corresponding to pressure range from 1 to 1.5 x 10 -5 
atm. We have found that the ratio V/5 is independent of the altitude and then of pressure 
in the range studied. 

Moreover we have performed simulations to study the influence of the initial seed sep- 
aration yo on T>/5. Figure [5] shows the dependence of V/5 for yo < Q.2 No/N cm. Two 
sets of seed parameters were considered: (a) for yo > 0.1 Nq/N cm: n max = 10 13 N 2 /Nq cm, 
a = 0.02 N /N cm, and E = 1.5£ bd ; (b) for y < 0.1 No/Ncm: n max = 10 W N 2 /N^cm, 
a = 0.00736 N /N cm, and E — 2.5i?bd- The first set corresponds to the reference condition 
used for Figure [TJ The second set was chosen to be close to the conditions used in [9] , and 
is characterized by a high value of the applied electric field, comparable to typical values of 
the electric field in streamer heads, and small initial seed width and separation. Because of 
the small seed width of the second set, for this set, computations were performed for a space 
resolution of 1.95 /im. Figure [6] shows that for positive streamers with yo > 0.04 No/Ncm, 
the value of the ratio V/5 for the second set is slightly less than the reference value of 0.35 
for the first set, which can be explained by the small dependence of V/5 on n max , as seen 
in Figure [5l For negative streamers, Figure [6] shows that for yo > 0.04 No/N cm, the ra- 
tio V/ 5 is independent of the value of yo- F° r initial separation of streamers smaller than 
yo < 0.04 No/ N cm the transition between repulsion and merging is occurring very close to 
the initial Gaussian seeds and therefore the region where p > p max /2 is indistinguishably 
joined with the charge lobes, typically observed near the initial germs (shown in Figure [2]). 
In this case, the streamer head appears to be significantly wider, and consequently the ratio 
V/5 increases as shown in Figure [6l It is interesting to note that in the two sets used for Fig- 
ure |6l very different values of a have been used. We have carried out additional simulations 
to verify that the value of a has a negligible influence on the ratio V/5. 

Finally, for simulations presented in this work, we have found that a stable value for 
the ratio V/5 is obtained for merging streamers for an initial separation yo which is smaller 
or comparable to the longest absorption length of photoionization for air (0.1408No/N cm). 
These results clearly show the significant role of photoionization on streamer merging. We 
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Figure 3: (color online) Time evolution of the ratio T>/8 for (a) two negative streamers with 
photoionization source term, (b) two negative streamers without photoionization source term, 
(c) two positive streamers with photoionization source term. Results for initial separation of 
the seeds yo = 0.15 No/N cm. and j/o = 0.2 No/N cm are shown. Other parameters are similar 
to Figure 1. Critical merging ratios are indicated with stripes. 
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Figure 4: (color online) Ratios V/S for positive and negative streamers for different values 
of the applied electric field E a / E'bd- Other parameters are similar to Figure [TJ 
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Figure 5: (color online) Ratios T>/8 for positive and negative streamers for different values 
of the maximum density n max of seeds. Other parameters are similar to Figure [TJ 
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Figure 6: (color online) Ratios T>/8 for positive and negative streamers for different initial seed 
distances y , with (a) n max = 10 13 N 2 /Nq cm, E = 1.5i?bd, cr = 0.02 N Q /N cm; (b) n max = 
10 10 N 2 /Nfi cm, E = 2.5E hd , a = 0.00736 N /N cm. Absorption lengths for photoionization 
at ground pressure are shown. 
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have carried out simulations for larger values of yo, but it was not possible to derive simple 
conclusions on streamer interactions from the analysis of the value of the ratio V/5. 

In conclusion, in this work, we have carried out a parametric study on the interaction 
of two streamers in air using 2D numerical simulations based on an accurate numerical 
method which ensures a time-space error control of the solution. This study is a first step 
towards a quantitative analysis of streamer merging in real 3D conditions. We have shown 
that for a initial separation between streamers that is smaller or comparable to the longest 
characteristic absorption length of photoionization (0.1408Nq/N cm), streamers will start to 
merge, i.e., the distance between their trajectories starts to decrease despite their electrostatic 
repulsion, at the moment when the ratio D/S of the streamer characteristic width V and 
the streamer mutual distance S attains a value of about 0.35 for positive streamers, and 
0.4 for negative ones. These ratios are independent of the applied field, the initial seed 
configuration, and pressure. Moreover presented results clearly illustrate the significant role 
of photoionization on the streamer merging. Further studies are required to find the key 
parameters to characterize streamer interactions when 3/0 is considerably larger than the 
longest photoionization absorption length. 
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